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ABSTRACT 



Aims. We present our analysis of a local {z = 0.04 ~ 0.2) sample of 31 galaxy clusters with the aim of measuring 

■ the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set 
constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. 

Methods. We have exploited the large field-of-view and low instrumental background of ROSAT /PSPC to trace the 
density of the intracluster gas out to the virial radius. We stacked the density profiles to detect a signal beyond r2oo and 
^ ■ measured the typical density and scatter in cluster outskirts. We also computed the azimuthal scatter of the profiles 

$_H ' with respect to the mean value to look for deviations from spherical symmetry. Finally, we compared our average density 

I and scatter profiles with the results of numerical simulations. 

, Results. As opposed to some recent Suzaku results, and confirming previous evidence from ROSAT and Chandra, we 

observe a steepening of the density profiles beyond ~ rsoo- Comparing our density profiles with simulations, we find that 
non-radiative runs predict density profiles that are too steep, whereas runs including additional physics and/or treating 
CO ' gas clumping agree better with the observed gas distribution. We report high-confidence detection of a systematic 

^ , difference between cool-core and non cool-core clusters beyond ~ 0.3r2oo, which we explain by a different distribution 

■ of the gas in the two classes. Beyond ~ rsoo, galaxy clusters deviate significantly from spherical symmetry, with only 

■ small differences between relaxed and disturbed systems. We find good agreement between the observed and predicted 
' scatter profiles, but only when the 1% densest clumps are filtered out in the ENZD simulations. 

, Conclusions. Comparing our results with numerical simulations, we find that non-radiative simulations fail to reproduce 

the gas distribution, even well outside cluster cores. Although their general behavior agrees more closely with the 
observations, simulations including cooling and star formation convert a large amount of gas into stars, which results in 
a low gas fraction with respect to the observations. Consequently, a detailed treatment of gas cooling, star formation, 
AGN feedback, and consideration of gas clumping is required to construct realistic models of the outer regions of 
clusters. 
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1. Introduction of the X-ray emitting gas in the outer regions of galaxy clus- 
ters is important for mapping the gas throughout the entire 

The outskirts of galaxy clusters are the regions where cluster volume, studying the formation processes currently 
the transition between the virialized gas of clusters ^O''^ ^ the Universe, and performing accurate inass es- 

and the accreting matter from large-scale structure oc- timates for cosmological purposes (e.g., [ Alien et al. || 201l D. 
curs and where the current activity of structure for- 
mation takes place. Around the virial radius, the as- 
sumption of hydrostatic equilibrium, which is a neces- 
sary assumption for reconstructing cluster masses from Because of the low surface brightness of the X-ray 

X-ray measuremen ts, might not be valid any more (e.g., emitting gas and the extended nature of the sources, 

lEvrard et al.l Il996[). which c ould introduce biases to X- measuring the state of the intracluster gas arou nd th e 



ray mass proxies faasia et al.l[200 4; Piffa rcttj fc Valdarnini virial radius is challenging (Ettori fc Molendil l201l[ ). 



2008; Nagai et al.ll2007b 



iLau et al.lr2009l : lMeneghetti et al 



Recently, the Suzaku satellite has achieved a break- 
2010^ .Fabian et al.ll201lh . As a result, the characterization through in this domain, performing measurements of 
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cluster temperatures out to r20(El (|ReiDrich et al. 20091 
Bautz et al. '2009'; 'Kawaharada et al."2010': 'Hoshino et al.l 



2010t ISi mioncscu et al. 2011; Akamatsu ct al. 20ill 



Humphrey e t al. 2012), even in one case beyond r2oo 
(|George et al., .2009i) . although thi s detection i s likel y 
hampered by systematic effects (| Eckert et al.l l2011aD . 
Interestingly, some of the Suzaku results indicate very 
steep temperature profiles and shallow density profiles 
in cluster outskirts, at variance with the results from 



XMM-Newton (iPratt et all 120071: iLeccardi fc Molendil 



120081: ISnowden et al.l 120081: ICroston et all 12008 
(IVikhlinin et al.l '2006'; 'Ettori fc Balestrai 1200^ 



( Vikhlinin 
(Vikhlinin et al 



2005) 



Chandra 
20091) ■ ROSAT 
and with the 



. 1999; NeumannI . _ 

results from numerical si mulations (IRoncarelli et al.l 120061 : 
iTozzi fc NormanI 120011: iNagai fc Laul |2011D . TbTs. the 
behavior of the gas in cluster outskirts is still the subject 
of debate. Throughout paper, we refer to cluster outskirts 
as the region with r > r^oo . 

Thanks to its large field of view (FOV, ^ 2 deg^) 
and low instrumental background, ROSAT /PSPC is to the 
present day the most sensitive instrument for low surface- 
brightness emission. Its ability to de tect cluster emis s ion at 
larg e radii has been de monstrated bv lVikhlinin et ahl (|1999[ ) 
and[N euman 3 (l2Q0l) (hereafter, V99 and N05). Because 
of the large FOV, it can perform simultaneous local back- 
ground measurements, so it is less affected than Suzaku by 
systematic uncertainties. Its main limitation, however, is 
the restricted band pass and poor spectral resolution, which 
makes it impossible to measure cluster temperatures. 

This paper presents the analysis of a sample of 31 
galaxy clusters observed with ROSAT /PSPC, with the aim 
of characterizing the cluster emission at large radii and 
comparing the r esults with three different sets of numer- 
ical simulati ons (I Roncarelli et al.ll2006l ; iNagai fc Laull20lit 
IVazza et al.l '2010). The paper is organized as follows. In 
Sect, m we describe our cluster sample and the available 
data. We present our data analysis technique in Sect. |3I and 
report our results in Sect. S) We compare our results with 
numerical simulations in Sect. [5] and discuss them in Sect. 

n 

Throughout the paper, we assume a ACDM cosmology 
with = 0.3, nA = 0.7, fib = 0.047, and Hq = 70 km s^^ 
Mpc"^ 



2. The sample 

We selected objects in the redshift range 0.04 — 0.2, such 
that ^200 is easily contained within the FOV of the instru- 
ment and is large enough to allow for an adequate sam- 
pling of the density profile. We restricted ourselves to ob- 
servations with enough statistics to constrain the emission 
around the virial radius. Our final sample comprises 31 clus- 
ters in the temperature range 2.5-9 keV, with the addition 
of A2163 {kT ~ 18 keV). Among our sample, we classified 
14 clusters as cool core (CC) following the classification 
of ICavagnolo et all ()2009l) (i.e. they exhibit a central en- 
tropy Kq < 30 keV cm^), and 17 as non cool core (NCC, 
Kq > 30 keV cm^). We recall that CC clusters exhibit a 
relaxed morphology, a high central density and a temper- 
ature decrement in the central regions, while NCCs trace 
dynamically-disturbed clusters with irregular morphologies 



and flat tempera,ture and density profiles in their cores (e.g., 
ISanderson et al.ll2009l: iHudson et al.ll2010t) . 

Our sample of clusters, together with the log of the 
available data and some important quantities, is shown in 
Table [31 In Fig. [T] we plot the distribution of temperature 
(left hand panel) and central entropy (right hand panel) 
for our sample. It should be noted that the sample was 
selected based on the quality of the existing observations 
and might be subject to selection effects. However, for the 
purpose of this work we did not require that the sample 
be representative or complete, since we are interested in 
characterizing cluster outskirts, which exhibit a high level 
of self-similarity. 



3. Data analysis 

3.1. Data reduction 

We used the ROSAT Extended Source Analysis Software 
(jSnowden et al.l IT994I ) for data reduction. We filtered out 
time periods when the master veto count rate exceeds 220 
cts/sec (using valid_times), and extracted light curves for 
the whole observation using rate_pspc. We used the ao 
executable to model the atmospheric column density for 
the scattering of solar X-rays, and fit the light curves in 
each energy band to get the relative contributions of the 
scattered solar X-rays (SSX) and of the long-term enhance- 
ments (LTE), using the rate Jit executable. 

We then extracted event images in each energy band 
and the corresponding effective exposure maps, taking vi- 
gnetting effects into account. We computed the contri- 
bution of the various background components, the LTE 
(lte_pspc), the particle background (cast_part), and the 
SSX (cast_ssx), and combined them to get a map of all 
the non-cosmic background components. 

3.2. Surface-brightness profiles 

The point-spread function (PSF) of ROSAT /PSPC 
strongly depends on angle, and ranges from ~ 15 arcsec on- 
axis to 2 arcmin in the outer parts of the FOV. Thus, the 
sensitivity of the instrument to point sources is higher on- 
axis, and a larger fraction of the cosmic X-ray background 
(CXB) is resolved. Consequently, when detecting sources in 
the image it is important to use a constant flux threshold, 
such that the same fraction of the CXB is resolved over 
the entire FOV and the value measured in the source-free 
regions can be used to subtract the background. We detect 
point sources using the program detect with a minimum 
count rate of 0.003 cts/sec in the R3-7 band (~ 3 x 10"" 
ergs cm~^ s~^ in the 0.5-2.0 keV band) to resolve the same 
fraction of the CXB over the FOV, and mask the corre- 
sponding areas. To compute surface-brightness profiles, we 
extract count profiles from the event images in the R3-7 
band (0.42-2.01 keV) with 30 arcsec bins centered on the 
surface-brightness peak, out to the radius of 50 arcmin. 
We divide each pixel by its corresponding exposure to ac- 
count for the vignet ting effects, following the procedure of 
lEckert et all (l2011bfFI . We perform the same operation for 
the background map and subtract the non-cosmic back- 
ground profile in each bin. 



^ We define as the radius witliin which A/(< rA)/-^'n'r% = 



http:/ /www. iasf-milano.inaf.it/~eckert/newsite/Proflit. html 
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Fig. 1. Distribution of temperature (left) and central entropy (right) of the members of our sample (see Table [3]). In the 
left panel, A2163 {kT ^ 18 keV) is located outside of the range. 



We tested this procedure on four different blank fields 
to estimate the accuracy in our determination of the CXB. 
We extracted the surface-brightness profile for the four ob- 
servations from the center of the FOV, grouped the bins 
to ensure a minimum of 100 counts per bin, and fitted 
the resulting profiles with a constant (see Fig. [5]). While 
the agreement is qualitatively good, significant deviations 
from the model are found, leading to an excess scatter of 
^ 6%, which we used as an estimate of the systematic un- 
certainties in measuring the CXB. This value encompasses 
both the cosmic variance and the true systematic uncertain- 
ties, e.g., in the vignetting correction or determination of 
the particle background. The higher level of scatter in the 
central regions is explained by the small area of the cor- 
responding annuli, which implies a large cosmic variance 
likely due to discrete sources with fluxes just under our ex- 
clusion threshold. Since in most cases the value of r2oo is 
larger than 15 arcmin, our systematic error of 6% is a con- 
servative estimate of the level of systematic uncertainties 
at the virial radius. 

For each cluster, we then use temperatur e pro- 
files f rom the lite rature (XMM-Ne w ton. . Snow d en et al.l 
(I2008D: Chandra, ICa vagnoloeLa]] (l2009|); BeppoSAX, 
iDe Grandi fc Molendil (120021 ) to estimate the virial tem- 
perature of the cluster. We approximated T^r as the mean 
temperature in the 200-500 kpc region, i.e. excluding the 
cool core and the tempe rature decline in the outskirts 
(jLeccardi fc Molendill2008[) . Using this estimate of T^r, we 
computed the value of r2oo from the scaling relations of 
lArnaud et al. (l2005l) . One migh t argue that the scaling re- 
lations of lArnaud et al.l (|2005D were computed using the 
mean temperature in the 0.1 — 0.5r2oo region, which in most 
cases extends beyond the available t emperature profiles. 
Using the mean temperature profiles of lLeccardi &: Molendil 
(|2008f ). we computed the mean temperature extracted in 
the 0.1 — 0.5r2oo and 200-500 kpc regions. In the temper- 
ature range of our sample, we found that the results differ 
at most by 2%, so our values of r2oo are unbiased. We then 
used the source-free region of the observation (r > 1.3r2oo) 
to fit the surface-brightness profile with a constant and get 
the cosmic background level for the observation, with the 
exception of the Triangulum Australis cluster, for which we 
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Fig. 2. Surface-brightness profiles for 4 blank-field PSPC 
observations from the center of the FOV, fitted with a con- 
stant. The dashed line shows the vignetting correction curve 
for comparison, in arbitrary units; the bump at ~ 22 ar- 
cmin is caused by the support structure. The bottom panel 
shows the ratio between data and model. 



used the range r > l.lr2oo because of the high value of r2oo 
37 arcmin). 

After having estimated the sky background for our ob- 
servation, we again extracted the surface-brightness profile 
in the radial range — 1.3r2oo with logarithmic bin size. 
The best-fit value for the CXB was subtracted from the 
profile and its error was added in quadrature to each bin. 
The systematic error of 6% on the CXB was also added in 
quadrature to account for the cosmic variance and system- 
atic uncertainties. For comparison, we note that in most 
cases the statistical uncertainties in the profiles are on the 
order of 10% of the CXB value around r2oo- 
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3.3. Density profiles 

To compute the density profiles, we first rebinned our 
background-subtracted surface-brightness profiles to ensure 
a minimum of 200 counts per bin and a detection signifi- 
cance of at least 3a, to reac h sufficient statisti cs in each bin. 
We used the procedure of iKriss et al.l ()1983D to deproject 
the observed profiles, and the PSPC response to convert 
the observed count rates into emission measu re, through 
the no rmalization of the MEKAL model (see lEckert et al.l 
l2011aL for details). 



Norm 



10 



-14 



A7:[dAil + z)Y 



(1) 



which is proportional to the emission measure. We assumed 
that the spectrum of our sources is described by an ab- 
sorbed MEKAL mod el with Nh fixed to the 2fcm value 
([Kalberla et al.l l2005[ ) and abundance fixed to 0.3^0. We 
used temperature profiles from the literature (see Table [3]) 
and interpolated them onto the same grid as the SB pro- 
files. The resulting model was then folded with the PSPC 
response, and the conversion from PSPC R3-7 count rate 
to emission measure was inferred. Beyond the limit of the 
temperature profiles, the temperature of the outermost an- 
nulus was used. We note that the conversion from PSPC 
count rate to emission measure is highly insensitive to the 
temperature; between 2 and 8 keV, the conversion factor 
changes at most by 4%. Once converted into the MEKAL 
normalization, we inferred the density profiles, assuming 
spherical symmetry and constant density into each shell. 

The error bars on the density profiles were estimated us- 
ing a Monte Carlo approach. In each case, we generated ID'* 
realizations of the surface-brightness profile using Poisson 
statistics, and performed the geometrical deprojection fol- 
lowing the method described above. The la error bars were 
then estimated by computing the root-mean square devia- 
tion (RMS) of our 10'' realizations of the density profile in 
each density bin. 

3.4. Azimuthal scatter profiles 

For the purpose of this work, we are also interested in the 
deviations in the X-ray emission from spherical symmetry. 
We divide our images into N azimuthal sectors with con- 
stant opening angle, and compute the surface-brightness 
profiles in each sector individually. We then compute the 
scatter of the various sectors with respect to th e mean 
profile, following the definition introduced by iVazza et al.l 

pTThh . 



£2 = 



1 ^ 

N ^ 



{SB, - {SB)y 
{SB)^ 



(2) 



where (SB) is the mean surface-brightness and SBi,i = 
1..N denotes the surface-brightness computed in the vari- 
ous sectors. It must be noted that the statistical fiuctua- 
tions of the SB between the different sectors introduce a 
certain level of scatter in Eq. [21 which must be taken into 
account for determining the level of intrinsic scatter. We 
used two different methods to disentangle between statis- 
tical and intrinsic scatter. In the first case, we computed 
the level of statistical scatter independently and subtracted 



it from Eq. [S] In the second case, we used a maximum- 
likelihood estimator to determine the intrinsic scatter and 
its uncertainties. The two methods gave consistent results 
and are described in detail in Appendix [21 For the remain- 
der of the paper, we refer to the results obtained using the 
direct method (see Sect. lA.l]) . 

In our analysis, we group the bins of the total surface- 
brightness profiles to reach a minimum of 8a per bin to en- 
sure adequate statistics in the scatter measurements, and 
then divide our images into 12 sectors with an opening of 
30°. The result of this analysis is a radial profile describing 
the intrinsic azimuthal scatter of the X-ray surface bright- 
ness, in percent. 

It must be noted that the method presented here is sen- 
sitive to all kinds of deviations from spherical symmetry, 
whether it is induced by the asymmetry of the large-scale 
structure (e.g., filaments), by gas clumping or by ellipticity. 
The cause of the observed asymmetry cannot be determined 
from the azimuthal scatter alone. 



4. Results 

4.1. Emission measure and density profiles 

In Fig.[31we show the scaled emission measure profiles (left, 
following Eq.[l]) and the deprojected density profiles (right) 
for the 31 clusters in our sample. A self-s imilar scaling 
was applied to the emission-measure profiles (jArnaud et al.l 
,2002 ): i.e., each profile was rescaled by the quantity 



9/2 



kT y/' 

10 kcV j 



(3) 



The density profiles were rescaled by E'^{z) = ftm{l + 
z)^ + fol lowing their expected evolution with redshift 
(ICrosto n et al. 2008). As a l ready noted by s everal authors 
(e.g., Vikhh ni'netrall 119991 : iNeumannI [20051: ICroston et aD 



l2008i : iLeccardi et al. 



201011 . the profiles show a remarkable 



level of self-similarity outside of the core (r > 0.2r2oo)- 
On the other hand, the large scatter observed in the cen- 
tral regions reflects the distinction of the cluster popu- 
lation into CCs, showing a prominent surface-brightness 
peak, and NCCs, which exhibit a flat surface brightness pro- 
flle in their cores, as expected fro m the standard /3-model 
(jCavaliere fc Fusco-Femianolll976[ ). 



SB{r) ^SBall 




-3P+0.5 



(4) 



In the radial range 0.2 — 0.7r2oo, the scatter of the den- 
sity profil es is 10%-20%, in exce llent agreement with the 
Chan dra (IVikhliniii etahl |2006[) an d XMM-Newton re- 
sults (jCroston et al.ll2008D . However, ICroston et all (|2008[) 
needed to rescale the profiles by T~^/^ to account for the 
lower gas fraction in low-mass objects. In our case, perform- 
ing such a scaling does not reduce the scatter of the profiles 
further. This is probably explained by the relatively narrow 
temperature range spanned in our sample (all but one ob- 
jects have a temperature higher than 3 keV), such that the 
clusters in our sample should show little dependence on gas 
fraction. 
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10' 



Fig. 3. Scaled emission measure (left, in units of cm ^ Mpc) and density profiles (right) for the 31 clusters of our sample 
(see Table O. 



4.2. Stacked emission-measure profiles 

To compute the mean profile of our sample, we interpolated 
each profile following a predefined binning in units of r2oo 
common to all clusters and performed a weighted mean to 
compute stacked profiles. The errors on the interpolated 
points were propagated to the stacked profiles. We also di- 
vided our sample into the two classes (CC and NCC) to 
look for differences between them. 

In Fig. |4] we show the stacked emission-measure (EM) 
profile for the entire sample compared to the profiles 
stacked for the two populations separately (see also 
Appendix [C]). Interestingly, we note a clear distinction be- 
tween the two classes in cluster outskirts (see the bottom 
panel of the figure). Namely, beyond ~ 0.3r2oo, NCC pro- 
files systematically ex ceed CCs. A similar e ffect has re- 
cently been noted by iMaughan et al.l ()201lD . who found 
a crossing of the average density profiles at a similar ra- 
di us, and also at a lower stat istical signifi c ance i n the works 
of lArnaud eTall (|2010[ ) and iPratt etall ()2010[ ). We stress 
that this effect is really a difference between the two classes; 
i.e. it is not introduced by a biased distribution of another 
quantity (such as temperature or redshift). Indeed, group- 
ing the profiles according to the temperature or the redshift 
did not show any particular behavior, which indicates that 
we are really finding an intrinsic difference between the CC 
and NCC classes. This result could follow from a different 
distribution of the gas in the two populations or from a 
higher clumping factor in disturbed objects (see Sect. [SI). 

Alternatively, the observed difference could be explained 
by an inaccurate determination of r2oo for NCC clu sters. 
Indeed, the scaling relations of lArnaud et al.l ([200^ were 
computed under the assumption of hydrostatic equilibrium, 
which is fulfilled better in CC clusters. This explanation 
is, however, unlikely. Indeed, to recover self-similarity, our 
value of r2oo should have been systematically underesti- 

1 /2 

mated by ^ 10% for NCCs, i.e. since r2oo T^,/^ the virial 
temperature of the NCC clusters should have been under- 
estimated by more than 20%. From mock Chan dra observa- 
tions of a sample of simulated galaxy clusters, iNaeai et al.l 
(|2007b[ ) have determined that the spectroscopic tempera- 
tures of unrelaxed clusters differs from that of relaxed clus- 
ters by ~ 5%, which is not enough to explain the observed 
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Fig. 4. Stacked emission measure profile (in units of cm~^ 
Mpc) for the entire sample (black) , and the two populations 
individually (CC, red; NCC, blue). See also Appendix [Cl 
The bottom panel shows the ratio between the CC and 
NCC populations. 



difference. It is therefore unlikely that such a large error on 
the virial temperature would be made. 

We fit the mean scaled emission-measure profiles from 
Fig. |4] with the standard /3- model (Eq. |4]), adding a second 
/3 component in the case of the CC clusters to take the 
cool core into account. The (double) /? model gives a good 
representation of the data in the radial range — 0.7r2oo 
(~ fboo), but significantly exceeds the observed profiles 
above thi s radius, in agreement w ith the results of V99, 
N05, and lEttori fc Balestral (|2009t ). For CC clusters, the 
best-fit model gives /3 = 0.717 ±0.005, while for NCC clus- 
ters we find /? = 0.677 ±0.002. Fitting the radial profiles in 
the range 0.65-1.3r2oo, we observe a significant steepening, 
with aslope 13 = 0.963±0.054for CCs and /3 = 0.822±0.029 
for NCCs. As explained above, the slope of the NCC profile 
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10-* 
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10' 



Table 1. Mean emission- measure and density profiles com- 
puted from our sample. 



Fig. 5. Average proton density profile for the entire sample. 
The dashed lines indicate the positive and negative scatter 
of the profiles around the mean value. 



is flatter than that of the CC profile beyond rsoo- The fits 
of the profiles in various radial ranges are reported in Table 
[2] to quantify the steepening. 

Given the limited number of objects in our sample, we 
have to verify that this result is not a chance realization. 
We fit all the emission-measure profiles at r > 0.3r2oo with 
a /3 profile, fixing the value of (3 to 0.7 and Vc to 0.12r2oo, 
and extracted the best-fit normalization for all profiles. We 
then sorted the normalization values into the CC and NCC 
classes, and performed a Kolmogorov-Smirnov test to de- 
termine the probability that they originate in the same par- 
ent distribution. Using this procedure, we found that the 
chance probability for this result is very low, P^6xl0^^. 
Therefore, we can conclude with good confidence that we 
are indeed finding an intrinsic difference between the two 
classes. 
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0.93 


1.05 


(1.06 zb 0.11) 
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(6.33 ± 1.01) 
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22 



Note: Column description. 1 and 2: Inner and outer bin radii in 
units of r-2oo; 3: Emission measure rescaled by A^sc in units of 
cm~® Mpc; 4: Average proton density in units of 10"'^ cm"''; 
5: Scatter of the various profiles relative to the mean value in 
percent. 

Our density profiles are in good agreement with the re- 
sults of V99. However, while V99 estimated the density 
from /3-model fitting, we performed a geometrical depro- 
jection of the data using temperature profiles to infer the 
mean density profile. This method has the advantage of not 
depending on any model. 



4.3. Stacked density profiles 

We stacked the density profiles shown in the right hand 
panel of Fig. [3] following the same method as for the EM 
profiles. From the different profiles, we also computed the 
scatter of the profiles around the mean value, following a 
method similar to the one presented in Sect. 13.41 for the 
azimuthal scatter. The statistical scatter was subtracted 
from the total scatter using the same technique. In Fig. 
[5] we show the average density profile of our clusters to- 
gether with the scatter of the individual profiles around 
the mean value (see also Tabic [1]). At r2oo, the mean den- 
sity is 71200 = (3.8 ± 0.4) X 10-^£;2(z) cm'^, with 25% 
scatter. For comparison, it is interesting to note that the 
density of PKS 074 5-191 claimed in the Suzaku analysis of 
iGeorge et al.l (j2009l) at r2oo deviates from our mean value 
by more t han 5a, which casts e ven more doubt on this mea- 
surement (jEckert et al.ll2011aD . 

As for the EM, we also extracted mean density profiles 
individually for the two classes of clusters in our sample. 
The same behavior is observed at large radii; i.e., the den- 
sity of NCC clusters is systematically higher (by ~15%) 
than that of CCs above r ~ 0.3r2oo- A global steepening of 
the density profiles is also observed beyond ^ ^500- 



4.4. Gas mass 

We computed the gas mass from our deprojected density 
profiles and stacked them in the same way as described 
above. In the self-similar model, the gas mass is expected 
to fo llow the relation M oc T'^/^ (e.g., iBrvan fc Normaiil 
Il998l ). However, observational works indicate that the ac- 
tual Mgas — T relation is steeper than the expected self- 
similar scaling ([N euman n fc ArnaudI I2OOII : lArnaud et al.l 
120071 : ICroston et al.ll2008t) because of the lower gas fraction 
in groups and poor clusters. For this work, we use the rela- 
tion d etermined from the REXCESS sample (jCroston et al.l 
I2OO8D to rescale our gas mass profiles, 

"--^'''"(toIJv) ■ 

As above, we divided the sample into CC and NCC 
classes, and stacked the two classes individually. In Fig. [Hj 
we show the mean gas mass profiles for CC and NCC clus- 
ters. As expected, CCs have a higher gas mass in their inner 
regions, since their central densities are higher. More inter- 
estingly, we see that the two profiles converge in cluster 
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Fig. 6. Enclosed gas mass profiles for CC (red) and NCC 
systems (blue). The data were resca led by Eiz)kT~^-^^^ 
observed in the REXCESS sample (jCroston et al.ll2008[) . 



as 



outskirts, and exhibit a gas mass around the virial radius 
that is consistent within the error bars. At r2oo, the univer- 
sal gas mass is 



M, 



gas, 200 



(2.41±0.05)xl0''*£;(z) 



kT 



10 keV 



1.986 



Mo, (6) 



with a scatter of 17% around the mean value. This result 
follows from the higher density measured in average beyond 
~ 0.3r2oo in NCC clusters and the steeper slope of CC pro- 
files in the outskirts (see Sect. 14. 2p . The lower density of 
CC clusters in the outer regions compensates for the well- 
known excess observed in the cores, such that the total gas 
mass contained within the dark-matter halo follows a uni- 
versal relation. We also estimated the average gas fraction 
by computi ng the expecte d value of M200 using the scaling 
relations of lArnaud et al.l (|2005[ ). For our sample, we find 
a mean gas fraction within r2oo of 



kT \ 

/,.,2oo = (0.15 ± 0.01) (^^j 



0.486 



(7) 



in good agr e ement with previous work s (e.g., 
IVikhlinin et all l200fil: iMcCarthv et~an 1200711 . which 
for the most massive o b jects corresponds t o ^ 89% of the 
cosmic baryon fraction (fjarosik eTaPIMll) . 

4.5. Azimuthal scatter 

Following the method described in Sect. 13.41 we computed 
the azimuthal scatter of the surface-brightness profiles for 
all the clusters in our sample, and rescaled the scatter pro- 
files by our estimated value of r2oo- We then stacked the 
profiles using the same procedure as described above and 
computed the mean azimuthal scatter. We recall that since 
the surface brightness depends on , the variations in den- 
sity are less important than the ones computed here. 

In Fig. [7] we plot the average scatter profile, compared 
to the mean value for CC and NCC clusters. The increase 
in the innermost bin is an artifact introduced by the small 
number of pixels in the center of the images, so it should 
be neglected. At small radii (r < 0.5r2oo) we find a clear 



Fig. 7. Stacked azimuthal scatter in surface-brightness for 
the entire cluster sample (black). The red and blue data 
represent the mean profile extracted from CC and NCC 
clusters, respectively. 



difference between CC and NCC clusters, which is eas- 
ily explained by the more disturbed morphology of the 
latter. In this radial range, CC profiles exhibit a scat- 
ter of 20-30%, which corresponds to density variations on 
the order of 10%, in g ood ag reement with the value pre- 
dicted bv lVazza et al.l ()2011bl ) from numerical simulations. 
Conversely, beyond r ~ rsoo, the profiles for CC and NCC 
clusters are similar, and indicate a high scatter value (60- 
80%). 

We investigated whether any systematic effect could af- 
fect our result in cluster outskirts, where the background 
is dominating with respect to the source. Indeed, in such 
conditions, the total scatter is dominated by the statistical 
scatter. In case the mean level of systematic uncertainties in 
the CXB reconstruction exceeds our adopted value of 6%, 
Ea. lA.2l immediatelv implies that the intrinsic scatter would 
be overestimated. The presence of both intrinsic and sta- 
tistical scatter could also introduce some covariance term, 
which is not taken into account in Eq. IA.21 To test this hy- 
pothesis, we ran a set of simulations including source and 
background, where we introduced a given level of intrinsic 
scatter for the source and a systematic error in addition 
to the Poisson statistics for the background. We then com- 
puted the intrinsic level of scatter following Eq. IA.2I Our 
simulations indicate that, even when increasing the level of 
systematic uncertainties to 12% of the CXB value, a signif- 
icant bias in the measured scatter only appears when the 
source-to-background ratio is close to the systematic uncer- 
tainties. Since, by construction, we never detect any signal 
when the source is less than '-^15% of the CXB value, our 
results are unaffected by these effects, and we can conclude 
with good confidence that the high level of scatter mea- 
sured beyond ^ t^qq is an intrinsic property of our cluster 
sample. 

In addition, we also tested whether the scatter for the 
two populations in the outermost regions could be affected 
by small-number statistics or driven by some particular ob- 
jects. Of the 31 objects in our sample, a measurement of the 
scatter at r2oo could be obtained for 23 of them (12 NCC 
and 11 CC). We used a jackknife method to test whether a 
single object dominates the results for any of the two pop- 
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ulations; i.e., we randomly exclude one or two profiles from 
the sample, recompute the mean profiles, and examine the 
distribution of the mean values. In both cases, the distri- 
bution of results is regular, which indicates that our results 
are not biased by a particular object. 

V99 also investigated the deviations from spherical sym- 
metry by measuring the value of j3 in six sectors in the ra- 
dial range r > O.Srigo, and concluded that the assumption 
of spherical symmetry is relatively well satisfied in cluster 
outskirts, at variance with our results (see Fig. [7]). However, 
when fitting a /3-model the fit is mostly driven by the shape 
of the profile in the innermost region, where the statistics 
are higher. Conversely, our method is model-independent, 
and it directly stacks the data at similar radii. For relaxed 
objects, our data also indicate little deviation from spheri- 
cal symmetry at r < rsoo, and a significant scatter is only 
observed beyond rsoo, so it is probable that these devi- 
ations would not be reflected in the /3-model fit. For in- 
stance, the case of A2029 is striking. While, in agreement 
with V99, we find little azimuthal variations in Pouter, we 
observe a high level of scatter in this object beyond rsoo, 
which is explained by the presence of a possible filament 
connecting A2029 to its neighbor A2033 in the north (see 
iGastaldello et all (|2010[ ) and Appendix [B|). Moreover, V99 
deliberately excluded a number of systems with obviously 
disturbed morphologies, such as A3558 and A3266, which 
we included in our sample. Therefore, our results do not 
contradict those of V99. 



5. Comparison with numerical simulations 

In this section, we compare our observational re- 
sults with three differ e nt sets of numerica l simulations 
(iRoncarelli et all 120061: iNaeai et all l2007bl: IVazza et al.l 
I2OIOD . We analvze the results of a composite set of cosmo- 
logical runs, obtained by the different authors with slightly 
different cosmological and numerical setups. In addition, 
the preliminary data reduction was made on each dataset 
following independent post-processing techniques, aimed at 
assessing the role of gas clumping on the comparison be- 
tween simulated mock and real X-ray observations. Our 
aim in this project is to test the most general and converg- 
ing findings of such different runs against our observations 
with ROSAT /V^VG. 

5.1. Simulations 
5.1.1. ENZO 

We use a sample of 20 simulated clusters from the high- 
resolution and non-rad iative (NR) r e simul ations of massive 
systems presented in IVazza et al.l (|2010t ). In this set of 
simulations, adaptive mesh refinement in the ENZO 1.5 code 
(| Norman et al.l 120071 ) has been tailored to achieve high 
resolution in the innermost regions of clusters (following 
the increase in gas and DM overdensity) , and also in the 
outermost cluster regions, following the sharp fluctuations 
of the velocity field, associated with shocks and turbulent 
motions in the ICM. For a detailed presentation of the 
statistical properties of the thermal gas (and of turbulent 
motions) in t hese s imulat ed systems we refer the reader to 
IVazza et al.l (poTol [20TTal) . 



5.1.2. ART 

We analyze a sample of ten simulated cl usters with 
Tx > 2.5 keV from the sample presented in ' Naeai et al.l 
(|2007ai rbl). These simulations are performed using the adap- 
tive refinement tree (ART) N-body -|-gas-dynamics code 
(jKravtsovlll999l : iKravtsov et al.ll2002fl . which is a Eulerian 
code that uses adaptive refinement to achieve high spatial 
resolution (a few kpc) in self-consistent cosmological simu- 
lations. To assess the impact of cluster physics on the ICM 
properties, we compared two sets of clusters simulated with 
the same initial conditions but with different prescription of 
gas physics. In the first set, we performed hydrodynamical 
cluster simulations without gas cooling and star formation. 
We refer to this set of clusters as NR clusters. In the second 
set, we turned on the physics of galaxy formation, such as 
metallicity-dependent radiative cooling, star formation, su- 
pernova feedback, and a uniform UV background. We refer 
to this set of clusters as cooling-f-star formation (CSF) clus- 
ters. For detailed descriptions of th e gas physics and mo ck 
X-ray images we refer the r e ader t o iNagai et al.l ()2007al lbl). 

Following iNaeai fc Laul (|201lh . we also computed the 
dumping-corrected gas density profiles of X-ray emitting 
gas with T > 10^ K for comparisons with X-ray observa- 
tions. Indeed, the formation of dense clumps increases the 
emissivity of the gas, which leads to an overestimation of 
the measured gas density when the assumption of constant 
density in each shell is made. For these profiles, we com- 
puted the average squared density from the simulations in 
each radial bin and took the square root of the total to 
mi mic the reconstructio n of density profiles from real data 
fsee lNagai fc Laull201ll . for details). 

5.1.3. GADGET 

This set includes four mas sive halos sim ulated with the 
GADGET-2 Tree-SPH code (ISpringell l2005l). with M-.nn > 
lO^fM© (for a detailed description see IRoncarelli et al.l 
I2OO6I . and references therein). Each object was simulated 
following two different physical prescript ions: a NR run (re- 
ferred to as ovisc in Roncarelli et al.ll2006) and a run in- 
cluding cooling, star formation, and supernovae feedback 
(CSF). 

To eliminate the dense clumps that dominate the den- 
sity and surface brightness in the outskirts, when comput- 
ing the profiles for every radial bin, we excise the one per 
cent of the volume that corresponds to the densest SPH 
particles. This empirical method mimics the procedure of 
masking bright isolated regions from the analysis of ob- 
served clusters. 

5.2. Comparison of gas density profiles 

We compared the simulations with our observed mean 
ROSAT density profile (see Fig.[5]and TableP. We present 
the detailed comparison in Fig. [H with the NR simulations 
(left hand panel) and with the CSF simulations (right). 
From the figures, we find relatively good agreement be- 
tween all the different sets of simulations, especially beyond 
~ 0.7r2oo- The NR GADGET run has a lower normalization 
than the corresponding grid codes, because in GADGET the 
fraction of baryons virializing into clusters is less than the 
cosmic value (~ 78% of the cosmic baryon fraction), while 
grid codes predict a baryon fraction in clusters very close 
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Fig. 8. Comparison between the mean ROSAT density profile for our sample and the different sets of numerical simula- 
tions. The shaded area indicates the data and la scatter as shown in Fig. [Sj The bottom panels show the ratio between 
simulations and data as a function of radius. Left: Comparison with NR simulations . The dotted red cu rve represents the 
ENZO profile (Vazza et al. 2010), the solid green curve shows the ART simulations (jNagai et al.ll2007b[) . and the dashed 
blue curve is the GADGET profile (jRoncarelli et al.ll2006[ ). Right: Same with CSF simulations. The dashed blue line shows 
the GADGET simul ations, while the gr een curves show the ART profiles, for the total density (solid) and corrected for 
clumping (dotted. pNagai &: Laull2011l ). 



to the cosmic value. In general, we see that the predicted 
density profiles are too steep compared to the data. We 
note that NR runs predict steeper profiles than the runs in- 
cluding cooling, star formation, and feedback effects. CSF 
profiles also have lower normalizations, since radiative cool- 
ing transforms a fraction of the gas into stars. The profile 
including the effects of clumping shows the best agreement 
with the data. 

To quantify this effect, we fitted the various profiles in 
three different radial ranges (0.2 — 0.4r2oo, 0.4 — 0.65r2oo, 
and 0.65 — 1.2r2oo)- In the inner regions, the effects of addi- 
tional physics are expected to be important, thus highlight- 
ing the differences between NR and CSF runs. The radial 
range 0.4 — 0.65r2oo(~ 0.6 — Irsoo) is a good range for com- 
paring with the data, since the effects of radiative cooling 
should be small, and data from several different satellites 
are available for cross-check. On the observational side, the 
density profiles in this radial range are well-fitted by the (5- 
model (see Eq. 0]), and several independe nt works converge 
to the canonical value oi (3 ^ 0.7 (e.g., i Mohr et"allll999t 



Ettori fc Fabianlll999t IVikhhnin et al. 19991: ICroston et al.l 
20081: lEttori fc Balestral l2009t lEckert et all 1201 Ibfl . As a 

benchmark, we computed the values of /3 for our average 
density profile and the vario us sets of simulatio ns, fixing the 
core radius to 0.127-200 fe.g.. lMohr et al.iri999t) . The results 
of this analysis are shown in Table [2l The fits to the ob- 
servational data were performed on the emission-measure 
profiles (see Sect. 14. 2|) to take advantage of the larger num- 
ber of bins and minimize the uncertainties linked to the 
deprojection procedure. 

These numbers confirm the visual impression that the 
simulated gas density profiles are steeper than the observed 
ones. In the 0.4 — 0.65r2oo range, while all our datasets 
converge to a /3 value very close to the canonical value, 
all the simulations lead to significantly steeper gas profiles. 



Table 2. Values of the j3 parameter 

(jCavaliere &: Fusco-Femianol |1976[ ) in several radial 
ranges for the average ROSAT profiles and the various 
sets of simulations. 



Data set 


/3o.2-0.4 


/3o.4-0.65 


/3o.65-1.2 


Data, total 


0.661 ±0.002 


0.710 ± 0.009 


0.890 ± 0.026 


Data, CC 


0.700 ± 0.004 


0.699 ±0.016 


1.002 ± 0.057 


Data, NCC 


0.635 ± 0.003 


0.723 ±0.011 


0.853 ±0.029 


ENZO 


0.744 


0.945 


0.952 


ART, NR 


0.801 


0.956 


0.983 


ART, CSF 


0.808 


0.842 


1.005 


ART, NR, cl 


0.701 


0.824 


0.854 


ART, CSF, cl 


0.803 


0.718 


0.902 


GADGET, NR 


0.856 


0.857 


0.971 


GADGET, CSF 


0.756 


0.864 


0.944 



Note: The core radius was fixed to 0.12r2oo in all cases. The 
subscript cl indicates the profiles corre cted for the ef f ect of 
clumping using the method described in lNagai fc Laul (|2011l 'l. 



with (5 values higher than 0.85, with the exception of the 
ART profile that includes CSF and clumping. Therefore, we 
can see that at this level of precision the effects of additional 
physics cannot be neglected, even in regions well outside of 
the cluster core. 

The results presented in Table [5] also highlight the dif- 
ferences between NR and CSF runs. Inside rsoo, the sim- 
ulations including additional physics lead to flatter den- 
sity profiles compared to the NR runs. In this case, gas 
cooling converts a fraction of the X-ray emitting gas into 
stars. Since the cooling efficiency decreases with radius, 
more gas disappears from the X-ray range in the central 
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Fig. 9. Left: Comparison between the average observed azimuthal scatter profile from Fig. [7] (black) and the scatter in 
the simulations for the ENZO runs (red), for the total scatter (solid line) and when filtering out the 1% most-luminous 
cells (dashed curve). The cyan (NR) and magenta (CSF) curves represent the scatter in the ART simulations. Right: Same 
for the CC (red) and NCC (blue) observed profiles, compared to the 1%-filtered ENZO profiles for the morphologically 
relaxed (red) and disturbed (blue) simulated clusters. 



regions, which results in flatter density profiles and lower 
normalizations. We note, however, that this effect is prob- 
ably overestimated in the CSF simulations. Indeed, it is 
well-known that these simulations predict a stellar fraction 
that i s well above the observed value (e.g., Kravtso v et al.l 
l2005HBorgani fc Kravts"wll2009D . This effect is particularly 
strong in the ART CSF simulation, for which nearly one 
third of the gas is converted into stars. Beyond r^oo, there 
is little difference between NR and CSF runs; i.e., the effects 
of additional physics ar e not important. At large radii, the 
effect of gas clumping (jNaeai fc Laull20l"Th dominates and 
flattens the observed profiles. As we can see in Table [5] and 
in the right hand panel of Fig. [51 the ART profile including 
both additional physics and a post-processing treatment of 
clumping reproduces the behavior of the data more closely, 
even though it is still slightly too steep. 



5.3. Azimuthal scatter 

A study of the azimuthal scatter in the radial profiles of 
density, temperature, entropy and X-ray brigh tness of sim- 
ulated ENZO clusters has been presented in IVazza et al.l 
(l2011bD . In this case, we differ from the analysis reported 
there by computing the azimuthal scatter from more an- 
gu lar sectors, N=12, than for N=2, 4, and 8 explored 
in IVazza et al.l (|2011b[ ). In the simulations, several dense 
clumps are present, which may bias the predicted scatter. 
To overcome this problem, we computed the scatter of the 
simulated clusters both for the total gas distribution and 
by filtering out t he 1% most X-ray luminous cells, as in 
iRoncarelli et al.l ()2006l ). which removes a large fraction of 
the clumps. 

We also performed a similar analysis on the set of ART 
simulations, both for the NR and CSF runs. In this case, we 
analyzed mock X-ray images using the same method as the 
observational data (see Sect. 13. 4p . and applied our point- 



source detection algorithm to remove the most prominent 
clumps 

In Fig. [3J we show the measured scatter profile from 
Fig. [71 together with the scatter profiles of X-ray bright- 
ness from ENZO and ART simulations. Interestingly, we note 
that NR runs (red and cyan) overestimate the observed az- 
imuthal scatter, while CSF simulations underestimate it. 
In the latter case, radiative cooling lowers the entropy of 
the gas, which makes it sink into the cluster's potential 
well. This effect produces more spherical X-ray morpholo- 
gies, thus lowering the azimuthal scatter. Conversely, in NR 
runs, the effects of dynamics are more important, which cre- 
ate more substructures and increases the azimuthal scatter. 

Interestingly, the profile that best reproduces the data 
is the ENZO profile for which the 1% most-luminous pix- 
els were filtered out. This may indicate that some clumps 
are indeed present in the observations, but were detected 
as point sources and were masked for the analysis. We note 
that, even if in this case the azimuthal scatter from NR sim- 
ulation runs is in good agreement with the ROSAT data, 
the absolute profiles of density are too steep compared to 
observations (see the left hand panel of Fig. [5]). However, 
our definition of the azimuthal scatter (Eq. [2]) is normal- 
ized to the absolute value of the profile at each radii, which 
makes it a rather robust proxy of cluster asymmetries on 
large ~ Mpc scales. 

In the right hand panel of Fig. [HI we also show the aver- 
age radial trends of the azimuthal scatter for the projected 
X-ray emission from the ENZO clusters after dividing the 
dataset into 11 CC-like and 9 NCC-like objects, compared 
to the observed scatter profiles for the CC and NCC classes 
from Fig. [71 This division is of course only qualitative, since 
no radiative cooling is modeled in these runs. However, our 



Because of how few objects are considered, we ignored the 
GADGET simulations for this analysis. For a comparison be- 
tween GADGET and E MZD scatter profiles, we refer the reader to 
IVazza et al.l (|2011bD . 
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sample can be divided into classes that are quite similar 
to observed CC and NCC properties, based on the anal- 
ysis of the power ratios P3/P0 and of the centroid shift 
w, evaluated within rsoo as in ICassano et al.l (|2010l) . We 
classify as NCC-like systems those for which the values of 
P3/P0 > 10^^ and w > 0.02 were found in at least two of 
the three projected maps along the coordinate axes, or as 
CC-like o therwise, identical to what was done for the same 
sample in lVazza et al.l ()2011aD . 

In this figure, we can clearly see that the radial trend 
of the difference between the two populations disagrees. 
While in simulations the two trends detach as we move 
farther out in the ckister atmospheres, in the observed pro- 
files the most prominent differences are found in the range 
0.2 < r/r2oo < 0.8. In the CC case, we find better qualita- 
tive agreement in the outskirts than in the central regions. 
This is not surprising, given that radiative cooling and en- 
ergy feedback from central AGNs are missing in these runs. 
Indeed, as we can see in the left hand panel of Fig. [SI radia- 
tive coolin g has a stro ng impa ct on the general morphology 
of clusters ( Fang et al . 2009; Lau et ahl 2011). On the other 
hand, the simulated disturbed systems have a larger scatter 
in the outskirts than the observed NCC clusters. However, 
we observe large differences in the scatter between the var- 
ious NCC profiles, such that the result may be affected 
by small-number statistics. In any case, since the selection 
criteria are very different, we do not expect a one-to-one 
correlation between the various classes. 



6. Discussion 

6.1. Observational results 

In agreement w ith earlier works using R OSAT (V99, N05) 
and Chandra (|Ettori fc Balestral |2009[ ) , but at variance 
with some recent resul ts from Suzaku (iBautz et al.l l2009t 



ISimionescu et al.l I2OIII: iGeorge et al.l I2009D and XMM- 
Newton ([Urban et al.l 20 111) , our analvsis reveals that on 
average the slope of the density profiles steepens beyond 
'^500 (see Table [2|). This result indicates that the latter re- 
sults may have been performed along preferential directions 
connected with the large-scale structure (e.g., in the di- 
rection of filaments). Indeed, the narrow FOV of Suzaku 
only allowed sparse coverage of the outskirts of nearby 
clusters, so that these measurements might be the result 
of azim uthal variations. In the case of A1795. iBautz et aD 
([2009f ) detected a significant signal only in the northern di- 
rection, while the Perseus result (Simione scu et all 120111 ) 
was obtained along two narrow arms, covering less than 
10% of the cluster's extent at 7-200 • Moreover, using sev- 
eral offset ROSAT /PSPG pointings of the Perseus cluster, 
[Ett ori et al.l ([19981 ) observed clear azimuthal variations in 



the density and gas fraction. Therefore, it is likely that the 
aforementioned measurements are not representative of the 
cluster as a whole. This picture is supported by our analysis 
of azimuthal variations in cluster outskirts, which suggests 
that even CC clusters exhibit significant departures from 
spherical symmetry around r2oo- Consequently, a full az- 
imuthal coverage is required to study the global behavior 
of cluster outer regions. 

An important result of this work is the systematic dif- 
ference between CC and NCC cluster populations observed 
beyond ^ 0.3r2oo (see Fig. [J). As explained in Sect.[31 this 
effect seems to be an intrinsic difference between the two 



classes, since it is does not correspond to a biased distribu- 
tion of our sample in temperature or redshift. Our scaled 
gas mass profiles provide a natural explanation for this re- 
sult (see Fig. [5]). Indeed, when the appropriate scaling is 
applied, the steeper density profiles of CCs in the outskirts 
compensate exactly for the excess density in the central re- 
gions, such that clusters with the same virial mass have the 
same gas mass enclosed within r2oo, albeit distributed in 
a different way for relaxed and disturbed objects. This re- 
sult w as expected in the old cooling-flow scenario ([Fabianl 
I1994D . in which radiative cooling causes the gas to flow in- 
wards and accumulate in the central regions. While in the 
central regions AGN feedb ack prevents the gas from co oling 
below a certain level (e.g.. lMcNamara &: Nulsenll2007D . the 
entropy injected by the central AGN is not strong enough 
to balance the flow in the outer regions of clusters, which 
explains the steep density profiles seen in Fig.Uj Conversely, 
merging events are capable of injecting a very large amount 
of energy in the ICM, which results in an efficient redistri- 
bution of the gas between the core and the outer regions 
and creates the flatter density profiles measured for NCC 
clusters. 

We also determined the typical scatter in surface- 
brightness as a function of radius (see Fig. [7|) and split the 
data into the CC and NCC classes. In the central regions, 
we observe a systematic difference between CC and NCC 
clusters, with NCC clusters showing a higher level of scat- 
ter than CC. This result is easily explained by the larger 
num ber of substructures gen erally observed in NCC clusters 
fe.g.. [Sanderson et al. !'2009). For CC clusters, we measure 
a scatter of 20% — 30% below 0.5r2oo, which corresponds 
to small variations (~ 10%) in gas density. This indicates 
that the azimuthal scatter in the inner regions (r < 0.5r2oo) 
can be us ed to estimat e the X- ray state of clusters, as sug- 
gested by IVazza et al.l (|2011b( ) . Conversely, the scatter of 
CC profiles increases in cluster outskirts, and there is no 
observed difference between the two classes. Interestingly, 
we note that for CC clusters the turnover in Fig. [3 occurs 
around rsoo , which coincides with the radius beyond which 
large scale infall motions and filamentary accret ions are 
generally non-neghgible fe.g.. lEvrard et aD Il996[) . Inside 
^500 J the gas is virialized in the cluster's potential well, 
and it shows only little deviations from spherical symme- 
try. Beyond rsoo, accretion processes are important, and 
the gas is located mostly along preferential directions (i.e., 
filaments). As a result, the distribution of the gas becomes 
strongly anisotropic, even for clusters that exhibit a relaxed 
morphology in their inner regions. 

6.2. Comparison with simulations 

Comparing our density profiles with numerical simulations, 
we find that all NR simulations predict very steep profiles 
already starting from 0.2r2oo, with values of the /3 param- 
eter greater than 0.85 in the 0.4 — 0.65r2oo range (see the 
left hand panel of Fig. [8l and Table [2]) . This indicates that 
including non-gravitational effects is needed to reproduce 
the observed slope, even well outside of cluster cores. The 
runs including additional physics are in better qualitative 
agreement with the observations (see the right hand panel 
of Fig. [5]), although their gas fraction is too low because 
of overcoohng (^ 10% compared to ^ 15%). However, it 
seems unlikely that star formation and galactic winds (as 
in the CSF runs explored here) are the only feedback mech- 
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anisms needed to reproduce observed clusters. Indeed, sim- 
ple feedback models still face severe problems in matching 
the properties of the stellar components inside galaxy clus- 
ters, as well as the propert ies of galaxies within them (e.g., 
IBorgani fc Kravtsovll2009l for a recent review). 

As illustrated in Table HJ gas clumping may also play 
a role in reconciling simulations with observations. Indeed, 
if an important fraction of the gas in cluster outskirts is 
in the form of dense gas clumps, as suggested in simula- 
tions (jNagai fc Laull20ll] ). the emissivity of the gas would 
be significantly increased, thus leading to an overestima- 
tion of the gas density when the assumption of constant 
density in each shell is made. Our results show that the 
treatment of gas clumping slightly improves the agreement 
between data and simulations (see the right hand panel of 
Fig. [5]). In addition, gas clumping also provides an alterna- 
tive interpretation for our observed difference between the 
CC and NCC populations beyond 0.3r2oo- Indeed, simula- 
tions predict a larger clumping factor in unrelaxed clusters 
compared to relaxed systems for the same average density, 
which would result in a higher observed density in the for- 
mer. At the moment, it is not clear whether this difference 
is caused by gas redistribution or clumping, or if both of 
these effects play a role to some extent. 

On the other hand, we find that numerical simulations 
can reproduce qualitatively the observed azimuthal scatter 
in the galaxy cluster gas density profiles (see Fig. [S]), al- 
though they fail to reproduce the trends observed for the 
CC and NCC populations separately. Interestingly, we find 
that the observed azimuthal scatter is reproduced with rea- 
sonable accuracy when the 1% most luminous clumps are 
filtered out, whereas the NR simulations with no filtering 
overestimate the observed level of azimuthal scatter at all 
radii. Two possible interpretations can be put forward to 
interpret this result. Observationally, it is possible that the 
dense clumps were detected as point sources and were fil- 
tered out of our observations. If this is the case, long ex- 
posures with high- resolution X-ray telescopes ( Chandra or 
XMM-Newton) should allow us to characterize the point 
sources and distinguish between dense clumps and back- 
ground AGN, possibly unveiling the population of accreting 
clumps in cluster outskirts. Conversely, if such observations 
do not confirm the existence of the clumps, it would imply 
that NR simulations significantly overestimate the amount 
of clumping in cluster outskirts, which would weaken the 
case for the interpretation recently put forward to explain 
the flattening of the entropy profiles observ ed in a few cases 
(jSimionescu et al.ll201lUUrban et al.ll2011[ ). 

As shown in Fig. ^ radiative cooling may also help rec- 
oncile the NR simulations with the data. Indeed, radiative 
cooling lowers the entropy of the gas and makes it sink 
into the potential well, which produces clusters with more 
spherical morphologies ()Lau et al.ll2011[ ) and thus reduces 
the azimuthal scatter. Since we know that this effect is over- 
estimated in our CSF simulations, radiative cooling likely 
reduces the azimuthal scatter with respect to NR simula- 
tions, although not as much as what is predicted here. This 
effect may also explain why NR simulations fail to repro- 
duce the average scatter profiles of CC clusters (see the 
right hand panel of Fig. [S]). 

Alternatively, AGN feedback may be an important in- 
gredient that is rare ly taken into accou nt in numerical sim- 
ulations. Recently, iPratt etall (|2010[ ) observed an anti- 
correlation between entropy and gas fraction, such that 



multiplying cluster entropy profiles by the local gas frac- 
tion allows recovery of the entropy profiles predicted from 
adiabatic compression; i.e., the excess entropy observed in 
cluster cores is balanced by a lower gas fraction, and the to- 
tal entropy follows the p redictions of gravitational collapse. 
iMathews &: Gucl ()2011h interpret this result in terms of the 
total feedback energy injected in the ICM through various 
giant AGN outbursts, which they estimate to be as large 
as 10^'^ ergs. In this scenario, feedback mechanisms are pre- 
venting the gas from collapsing into the potential well, caus- 
ing a deficit of baryons in the inner regions of clusters, hence 
flattening the observed density profiles. Moreover, it is well 
known that this mechanism also takes place on group and 
galaxy scales, leading to shallower density profiles in the 
accreting clumps. As a result, the gas distribution in clus- 
ter outskirts would be more homogeneous than predicted in 
NR simulations, in agreement with our observed azimuthal 
scatter profiles. Therefore, although its implementation int o 
numerical simulations is challenging (jSiiacki et al.ll2008h . 
AGN feedback could be an important effect for reconciling 
simulations with observations. A more complex picture of 
the ICM, possibly including the detailed treatment of mag- 
netic fields, cosmic rays, and thermal conductions (and of 
the instabilities arising from these ingredients), would still 
represent a challenge for current cosmological simulations. 

7. Conclusion 

In this paper, we have presented our analysis of a sam- 
ple of local (z = 0.04 - 0.2) clusters with ROSAT /PSPC, 
focusing on the properties of the gas in cluster outskirts. 
We then compar ed our observational results with numeri- 
cal simulations ( Ronc arelli et al.l 120061 : iNagai fc Laul 120111 : 
Va zza e t al. 2011b). Our main results can be summarized 
as follows. 

— We observed a general trend of steepening in the ra- 
dial profiles of emission-measure and gas density be- 
yond 21 ''500 7 in good agre e ment with earlier works 

from Fvik hlinin et al. |1999D . iNeumannI (|2005D . and 
lEttori fc Balestra (,2009h . As a result, the shallow den- 
sity profiles observed in several clusters by Suzaku 
(iBautz et al.l l2009t ISimionescu et al.l l201ll ) are proba- 
bly induced by observations in preferential directions 
(e.g., filaments) and do not reflect the typical behavior 
of cluster outer regions. 

— We found that NCC clusters have in average a higher 
density than CC systems beyond ^ 0.3r2oo, which can- 
not be easily explained by any selection effect. We in- 
terpreted this result by a different distribution of the 
gas in the two populations: the well-known density ex- 
cess in the core of CC clusters is balanced by a slightly 
steeper proflle in the outskirts, which leads to the same 
gas mass enclosed within ^200 in the two populations 
(see Fig. . Alternatively, this result could be caused 
by a larger clumping factor in disturbed objects, leading 
to an overestimate of the gas density of NCC clusters 
in the external regions. 

— We also observed that NCC systems have higher az- 
imuthal scatter than CCs in the central regions, which 
is easily explained by the more disturbed morphology 
of NCC clusters. Conversely, beyond ~ rsoo both pop- 
ulations show a similar level of asymmetry (60-80%), 
which suggests that a signiflcant fraction of the gas is 
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in the form of accreting material from the large-scale 
structures. 

— Comparing our ROSAT density profile with numerical 
simulations, we found that all NR numerical simula- 
tions fail to reproduce the observed shape of the den- 
sity profile, predicting density profiles that are signif- 
icantly too steep compared to the data (see Table [5] 
and Fig. [SJ . This implies that nongravitational effects 
are important well outside the core region. The runs in- 
cluding additional physics (cooling, star formation, SN 
feedback) predict flatter profiles, although still too steep 
compared to the observations. Besides, it is well known 
that thes e simulations overpredict th e stellar fraction in 
clusters (jBorgani fc Kravtsovl |2009| ). A slightly better 
agreement is found when a treatment of the observa- 
tiona l effects of gas clumping is adopted (jNaeai fc Laul 

[mil). 

— NR simulations are able to predict the observed az- 
imuthal scatter profile with reasonable accuracy, but 
only when the 1% most luminous cells are filtered out 
(see Fig.[S]). This result implies that either (i) the clumps 
are quite bright and were masked as point sources in our 
analysis pipeline, in which case offset XMM-Newton and 
Chandra observations will be able to characterize them 
spatially and spectrally, or (ii) the non-radiative simu- 
lations significantly overestimate the effects of clumping 
on the observable X-ray properties. Because of the ab- 
sence of cooling, it is however hard for these simulations 
to reproduce the observed trends of azimuthal scatter 
for the two populations (CC and NCC) separately. 

As an alternative explanation, we suggest that AGN 
feedback might be important even at large radii, and could 
help to reco ncile observations and simulations. Indee d, re- 
cent works (Pratt et al.'' 2010l: iMathews fc Guoll201lD indi- 
cate that feedback mechanisms may be responsible for the 
well-known deficit of baryons in cluster cores, thus leading 
to flatter gas distributions out to large radii. Moreover, the 
existence of such mechanisms on group and galaxy scales 
could also dilute the accreting material at large radii, lead- 
ing to a smaller azimuthal scatter. 
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Table 3. Master table of the cluster sample. Column description: 1. Cluste r name; 2. Effective exposure of the PSPC observation; 3. Redshift (f rom NED); 4. 
Hydro gen column density, Nh, along the line of sight (jKalberla et al.l [20051) : 5. Mean temperature in the 200-500 kpc radial range; 6. rmn from lArnaud et all 



()2005f ) scaling relations, in physical units; 7. Same as 6, in appareii t un its; 8. Central density n p (t his work); 9. Central entropy -RTp, from lCa.vagnolo et al. ()2009[ ): 



. _ . _ W , 

10. Reference for the temperature profile (l= ISnowden et all (|2008^ : 2= ICavagnolo et al.l (|2009l ): 3= lDe Grandi fc MolendU (|2002D ) 
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Appendix A: Determination of azimuthal scatter 
profiles 

The azimuthal scatter (jVazza et alJ l2011lj ) is defined as 
the relative scatter in surface brightness between various 
sectors (see Sect. I3.4p . 



£2 = 



1 ^ 
N ^ 



{SB, - {SB)f 
{SB) 



(A.l) 



In practice, computing this quantity is difficult, since the 
statistical fluctuations of the surface brightness introduce a 
contribution to the scatter that is actually dominant in the 
outer regions. To estimate the intrinsic level of azimuthal 
scatter, we used two different complementary methods, 
which we describe in more detail here. 

A.l. Subtraction of the statistical scatter 

Since the statistical fluctuations of the data also introduce 
a certain level of scatter, it must be noted that the quantity 
computed through Eq. lA.ll gives the sum of the statistical 
and intrinsic scatter. 








6 


3 


< 


1 



1.2 



Fig. A.l. Comparison between the mean azimuthal scat- 
ter proflles computed using the direct method (black, see 
Sect. lA.ip and the alternative method using a maximum 
flkelihood estimator (red, see Sect. IA.2)) . 



The statistical scatter T,stat is given by the mean of the 
individual relative errors, 

1 ^ ^2 

SLt ^ J^Yl (^SB)'^' '•^'^^ 

and must be subtracted from Eq. [5] to estimate the level of 
intrinsic scatter. The validity of Eg. IA.3l for the statistical 
scatter was verified through a set of simulations of a source 
with no intrinsic scatter. 

The uncertainties in the scatter are then estimated 
through Monte Carlo simulations. Namely, the surface- 
brightness values in the N sectors are randomized, and the 
scatter is recomputed each time. This procedure is applied 
10'^ times, and the error on the scatter is defined as the 
RMS of the distribution around the mean value. 

A. 2. Maximum likelihood estimation 

To check the vafidity of our approach we performed an inde- 
pendent analysis of the scatter. We model the intrinsic scat- 
ter in the f orm of a Gaussian. We use a maximum likelihood 
algorithm ()Maccacaro et al.]|1988l) to fit the data, where the 
free parameters are the mean and the intrinsic scatter (i.e. 
the standard deviation of the Gaussian). The methods de- 
scribed in both Sect. lA.ll and this appendix were applied 
to the surface brightness distribution within the annuli of 
each cluster (see Sect 13.41 for details) . Intrinsic scatter pro- 
files from different objects were rebinned onto a common 
grid in units of r2oo and stacked. In Fig. lA.ll we compare 
the intrinsic scatters measured with the two methods. The 
profiles are very similar, the general trend towards increas- 
ing scatter with radius is recovered with both methods. The 
only bin where a significant difference is observed is around 
0.7r2oo. This comparison therefore provides a confirmation 
of our scatter analysis using two very different methods. 



Appendix B: Notes on individual objects 

- A85: 

A subcluster located ^ 10' south of the cluster center is 
currently merging with the main cluster. This substruc- 
ture was masked for the analysis. 
A401: 

The cluster is connected through a filament to its neigh- 
bor A399, located ^ 35' south-west of the center of 
A401. We extracted the surface-brightness profile in a 
sector of position angle 340-250° to avoid any contami- 
nation of A399 to our measurement of the CXB. 

- A478: 

The combination of a favorable temperature/redshift 
and a good-quality ROSAT observation allows us to 
reach the highest signal-to-noise ratio in the sample at 
r2oo for this strong CC cluster. As a result, the data from 
this cluster may contribute strongly when a weighted 
mean is performed. 
A6U: 

This NCC cluster exhibits an unusual decreasing az- 
imuthal scatter profile, showing large (close to 100%) 
scatter in its central regions, but no significant scatter 
around r2oo- 

- A2029: 

A probable filament connects A2029 to A2033, located 
~ 35' north of the center of A2029. The surface- 
brightness profile was extracted in a sector with position 
angle 140-80° to measure the CXB level. 

- A2142: 

Several PSPC observations of this famous cold-front 
cluster exist. For this work, we used the longest avail- 
able observation, which was pointed 16' south of the 
center of A2142. This is the only case in the sample for 
which the observation was not pointed on the target. 
A3 558 and A3562: 

These two clusters are located in the Shapley super- 
cluster and connected by a filament. Consequently, they 
show an unusually high azimuthal scatter in the out- 
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skirts. The CXB level was estimated by excluding the 

direction of the filament. 

A3667: 

This very disturbed cluster shows the highest emission- 
measure and density in the sample beyond ^ 0.2r2oo, 
and hence it could bias our average profiles, in particu- 
lar when computing the difference between the CC and 
NCC classes. However, removing it from the sample did 
not lead to any significant difference, either quantitative 
or qualitative. 
A4059: 

This is the most azimuthally-symmetric cluster in the 
sample. The azimuthal scatter for this cluster is consis- 
tent with at all radii. 
— Hydra A: 

A tail of emission (filament?) extends out to ^ 20' 
south-east of the cluster core. This leads to a very high 
azimuthal scatter (> 100%) around r2oo- 



Appendix C: Mean emission-measure profiles 

In Table [UTT] we give the mean self-similar scaled emission- 
measure profiles for the CC and NCC classes and the whole 
sample, as shown in Fig. 31 



Appendix D: Computing the gas fraction from 
density profiles 

The gas fraction in the observations and in the simulated 
clusters within an overdensity A can be computed directly 
from the profiles presented in Fig. [5] Indeed, by definition. 

Ma = Ap^^t^nrl, (D.l) 
where pcrit = |^ = 9.2 x 10'^° g cm'^. Then, 

Making the substitution x = we find the convenient 
formula 

3 

/ffas.A = -r / Pgas{x)x'^ dx. (D.3) 

^Pcrit Jo 



17 



D. Eckert et al.: The gas distribution in the outer regions of galaxy clusters 

Table C.l. Data of Fig. H) mean self-similar scaled emission-measure profiles for the whole sample and for the CC and 

NCC classes, in miits of cm~^ Mpc 
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